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Some Cases of Wave Motion due to a Submerged Obstacle. 

By T. H. Havelock, F.RS. 

(Eeceived May 14, 1917.) 

1. As far as I am aware, only one case of wave motion caused by a 
submerged obstacle has been worked out in any detail, namely the two- 
dimensional motion due to a circular cylinder ; for this case, Prof. Lamb has 
given a solution applicable when the cylinder is of small radius and is at a 
considerable depth.* The method can be extended to bodies of different 
shape, and my object in this paper is to work out the simplest three- 
dimensional case, the motion of a submerged sphere. 

The problem I have considered specially is the wave resistance of the 
submerged body. In the two-dimensional case, this is calculated by considera- 
tions of energy and work applied to the train of regular waves. But for a 
moving sphere the wave system is more complicated, like the well-known 
wave pattern for a moving point disturbance, and similar methods are not so 
easily applied ; I have therefore calculated directly the horizontal resultant 
of the fluid pressure on the sphere. Before working out this case, the 
analysis for the circular cylinder is repeated, because it is necessary to carry 
the approximation a stage further than in Prof. Lamb's solution in order to 
verify that the resultant horizontal pressure on the cylinder is the same as 
the wave resistance obtained By the method of energy. 

The stages in approximating to the velocity potential may be described in 
terms of successive images ; the first stage fa is the image of a uniform 
stream in the submerged body, the second stage fa is the image of fa in the 
free surface, the third fa is the image of fa in the submerged body, and so 
on. In order to keep the integrals convergent, a small Motional coefficient is 
introduced in the usual manner ; after the calculations have been carried out, 
the coefficient is made zero. Further, the solution for uniform motion is 
built up so that expressions can be found for the velocity potential at any 
time after the starting of the motion, although only the final steady state has 
been studied in detail. The wave resistance of a sphere is found to have the 
form const, x a^ 2 e~~ a ^ 2 Wi, i(a), in which a is 2gf/c 2 , with /the depth of the 
sphere and c its velocity ; Wi, i (a) is a confluent hypergeometric function. 
In order to graph the wave resistance as a function of the velocity, 
expansions have been found for this particular variety of the function 

* H. Lamb, ' Ann. di Matematica,' vol. 21, p. 237 ; also ' Hydrodynamics,' 4th ed., .p. 401. 
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"W*, m (a) ; it belongs to the logarithmic case for which a general expansion is 

not available. 

In general form the graph of the resistance is very similar to that of the 

circular cylinder. 

Circular Cylinder. 

2. The steady state for uniform motion of the cylinder may be attacked 
directly, as in Prof. Lamb's solution, but we shall adopt his suggestion of 
building it up from simple oscillations. : Take the axis of x in the free 
surface of the water, and the axis of y vertically upwards. A circular 
cylinder, of radius a, is making small oscillations parallel to Ox with velocity 
c cos at, the axis of the cylinder being horizontal and perpendicular to Ox, 
and the mean position of the centre being the point (0, — /). A first 
appproximation when the depth / is sufficiently large is found by ignoring 
the surface effect altogether and putting 

. <£ = ca 2 (% J r 2 ) e iat ; r 2 = x 2j t{y J tf) 2 . (1) 

This satisfies the boundary condition at the surface of the cylinder. For the 
next step, add a term Xi to the velocity potential so as to satisfy the 
conditions at the free surface, but ignoring meantime the disturbance 
produced thereby at the surface of the cylinder. The term Xi must be a 
potential function and it must satisfy the condition for deep water, namely, 
dX.i/dy = f or y = -co ; these conditions are fulfilled by 

J QO 
a (re) e«y sin kx die, (2) 

where a, is a function of k to be determined. This form is chosen because we 
can satisfy the conditions at the free surface by using an equivalent form 
for (1), since 

foo 
e-*to+f) sin Kxdic ; y+f> 0. (3) 

o 

The surface elevation is expressed similarly by 



1} = e i(rt 



CO 

/3 (k) sin kx die. (4) 



o 



In order to keep the various integrals convergent, we assume that the 
liquid has a slight amount of friction proportional to velocity ; in the sequel 
the results are simplified by making the frictional coefficient & tend to zero. 
In these circumstances the pressure equation is 

■p/p = const, —gy + ^-iq 2 . (5) 

Hence the conditions at the free surface are, neglecting the square of the 

velocity, 

■—gy + ^fy = const. ; —dcji/dy = dr)/dt. 
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Here $ is the velocity potential after (2) has been added to (1) ; thus the 
equations for a and f3 are 

ca 2 /ce- Kf —fcu = icr/3 
From these we obtain the expressions for Xi and tj, namely 



V (6) 



-Xi = ca 2 e i(rt 



o gK — o-^ + lfMr 



r] = ca 2 e i(xt ~ 2 ^-y^-2 e~ K f sin fcxd/c. (8) 



o g/c — a 2 -\-i/jLcr 
The expression for Xi can be divided into two parts 



Xi = — ca 2 e ifrt 



if ,a • 7 o o „• * f 00 ^"^-^sin ##efo: 



2 ■ • (9) 

Jo O — %fJbG—gK 



If we regard X x as the image of the oscillating cylinder in the free 
surface, we see from the form of the first integral in (9) that part of the 
image is a negative doublet at the image point (0,/). We obtain next the 
velocity potential of the motion produced by a sudden small displacement of 
the cylinder, and we take this to be equivalent to a momentary doublet of 
constant strength. Suppose then that at a time r a doublet is suddenly 
created, maintained constant for a time St, and , then annihilated. The 
velocity potential at any subsequent time t is given by a Fourier synthesis of 
the preceding results for an oscillating cylinder, and we have 



S 



T 



roo 



</> = — e^-^Wda-y (10) 

7T Jo 

where [(f)] is the sum of (1) and (9), omitting the factor e iot . 

Carrying out this integration for the value of cj> in (1) and for the first part 
of (9) gives simply the momentary doublet at the centre of the cylinder and 
the negative doublet at the image point. These doublets last for a short 
time St ; the subsequent fluid motion is contributed by the second part of 
(9). For this we have to evaluate the real part of 

— . da: t-~T>0. (11) 

o o^ — ipcr—gic 



\ 



We obtain the value by contour integration ; further we simplify the 
result by neglecting ju, 2 . We shall make /jl zero ultimately, but we must 
retain it sufficiently to keep the integrals convergent ; however, at one or 
two stages, superfluous terms may be omitted when it is clear that the final 
limiting values will not be affected. We find for (11) the value 

— ir/cYe-^V-^sm {?cY(t — r)}, 
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writing V for \/(g/fc) whenever it serves to simplify the notation. Hence 
the velocity potential of the subsequent fluid motion after the cylinder has 
been given a small displacement at time r is 

= 2ca 2 Sre~^ t ~^ fcYe~*V-rt sin kx sin *V (t-r) die. (12) 

Finally we obtain the velocity potential for a cylinder in uniform motion 
by substituting x-\-c(t — r) for x, noting that hereafter x will refer to a 
moving origin immediately over the centre of the cylinder; we then 
integrate with respect to t from the start of the motion up to the instant in 
question. We could in this way obtain results for any stage of the motion, 
but we limit the discussion to the final steady state ; for this we take — cc as 
the lower limit in integrating with respect to r. Before writing down the 
result, we must remember to introduce the integrated effect of the original 
momentary doublet in (1) and its negative image, which were not included in 
(11) ; these clearly add up to steady doublets. Hence we find for the steady 
state 

poo 

<p = D — D 1 + 2ca 2 e- K (f~ri(Asm/cx + RcoSfcx)dic, (13) 

where D represents the doublet ca?x/r 2 at the point (0, — /), Di an equal 
doublet at the point (0,/), and 



**V(V + c) , /c 2 V(V-~c) 



2A = - v vv ^ v; + 

O /TT , \0 - 1 O ' 



2 J 



tc 2 (Y + cy + i^ « a (V-c)* + i,* 

3. Before proceeding further we may obtain the surface elevation from (13) 
for comparison. The surface condition is now 

—dfyjdy = dr)/dt = —cdrj/cix. 
Hence we have 

Too- 

7j = 2a 2 f/(x 2 -\-f 2 )—2a 2 /e (A cos /ex— B sin /ex) e~ K fdtc, (15) 

in which /e = g/c 2 . Further since yu, is to be small, we may omit irrelevant 
terms and put 

A = —Ko(/e--/eo)/{fe--(ieQ-{-ifj,/c)}{/e--(ieo--ifjL/c)}, 

B = fCo(/J>/c)/{/e — (/co + iji/c)} {k — (k — ip/c)}. (16) 

The integral in (15) can then be written as 

J ! +■ e Ve-'tdK. (17) 

\^K — Ko — t/MjC K — Ko + lfA/Cj 
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We transform these integrals by contour integration in the plane of a 
complex variable k } treating separately the cases of x positive and x negative ; 
after making ^ zero in the final results we obtain 

2a 2 f „ , . f°°mcosm/— /cosinm/ w ^ 7 . 

xr+f Jo m J + /t J 

2a 2 / rt o f 00 m cos m/—r sin mf 



V- 9 — jL e- mx dm: x>0. (18) 

mM-# 7 



These agree with Lamb's results for the circular cylinder in a uniform stream. 

The wave resistance E is derived from the regular waves in the rear, by 
considering the rate of increase of energy and taking into account the 
propagation of energy in a regular train ; we have 

K = \gp (amplitude) 2 = 4,ir 2 gpa 4: tc^e"~ 2,i ^. (19) 

4. We have now to obtain the resistance R by direct summation of the 
horizontal component of fluid pressure on the cylinder. It is clearly 
necessary to proceed to a further stage with the velocity potential, since we 
have assumed so far that the surface effect is negligible in the neighbourhood 
of the cylinder. If we write (13) as 

<I> = T> + X 1 , (20) 

the doublet D is the first approximation, satisfying the boundary conditions 
on the cylinder ; Xi is' the image of the doublet in the free surface, found by 
satisfying the conditions there. The next step is to find X 2 , the image of Xi 
in the cylinder, ignoring then the effect of X 2 at the free surface. It follows 
that X 2 is the image of Xi in the cylinder, found as if the cylinder were at 
rest in a field defined by X x . Taking polar co-ordinates with the origin at the 
centre of the circular section of the cylinder, we have 

x~rcos0; y+f = rsin.0 ; (21) 

also the conditions for X 2 are that it should be a potential function, the 
components of velocity must vanish as r becomes infinite, and 

3 (Xx + X a )/3r = 0, for r = a. (22) 

But from (13), Xi consists of a summation of terms of the fori. 

.,,,, COS 

e^y . /cos. 

sm 

We obtain X 2 by replacing each term by the expressions 

6 -~ K f eK a*smO/r COS , 2 ^ q'., 

sm v ' 7 

and the above conditions for X 2 are then satisfied. This process amounts 
simply to inversion : we may think of Xi as due to a line distribution of 
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sources and X 2 is then a circle of sources on the inverse of this line with 
respect to the cylinder. We have now for the velocity potential to this stage 

e -x if-y) { (A — | ) sin /ex + B cos kx) dtc 





o 

'00 



+ 2w 







e-Kf+^jM { (A- J) sin (fca 2 xjr 2 ) + B cos (tca 2 x/r 2 ) } d/c. (23) 



We have put A — J for A so as to include under the integral sign the 
doublet previously denoted by Di. 

The method could theoretically be carried on step by step ; however, we 
stop at this stage because it is sufficient for obtaining the wave resistance 
E from the pressure equation to the same approximation as by the energy 
method. 



We have E = 



c2tt 



ap cos dO ; (24) 



p/p = —cdcj>/dx—gtf + ^—^q 2 . (25) 

If we write (23) as <£ = D + Xi + X 2 , and omit terms which obviously 
contribute nothing to the value of E, we have, when r = a, 

p ex a 2 30 00 

= (2c/a)sinfl3(Xi + X 2 )/a« + ^(X 1 + X a ;, (26) 

where we have used (22) and the value of D. From (23), omitting the 
doublets D and Di, which will from symmetry give no contribution to E when 
jjl is zero, we have 



fXt 



p = 4ea 2 e~ 2 «f + ™ ™i o 1 2/cc A sin 6 sin (<£ - 0) + fi A sin $ 

+ 2/ccB $in#cos(<£ — 0) + /nBGO^^>}d/c i (27) 

where cj> = /ca cos 0. Substituting in (24) we have an expression for E. We 
may now change the order of integration and take first that with respect 
to ; we can carry this out, after some transformation, by means of the 
integrals 

* e h cos e cos (h sin — n0) d0 = 7Th n /T(n + l) } 



e h cos e cos ty s i n q + n Q^ de __ 0? ^28) 

where n is a positive, odd integer. In fact the integration with respect to 
gives simply 7r/ta(fccB+fiA)', hence we have 



'00 



E = 4fwpccfi K(fccB + jnA)e~ 2fi fdfc, (29) 
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where A and B are given by (14), or by (16) since we suppose /n small. Thus 
we have 

li = 4z7roca* Lim , 7 — •— ; — , x -. , 7 , , v ^^ 

^_>oJo {^ — (*o + W c )H*~- (/c — W c )} 

= 4zTTpccfi Lim yu- {2m/co 3 e~ 2K( > f / 2i {[ije) + finite quantity} 
= 4>n*tfpa*c-*e-WI* 9 (30) 

which is the same as the previous expression (19). 

Sphere. 

5. A sphere of radius a is at depth / below the surface and is moving with 
uniform velocity c parallel to the axis of x. The origin is in the free surface, 
the axis of z being drawn vertically upwards. As before, the first approxi- 
mation is a doublet D given by 

c/> = ca?xj 2r 3 ; r 2 = x 2 + f + (z +/) 2 . (31) 

For the purpose of satisfying the conditions at the free surface we have 

= D=-Jca 3 ~ e- K ^ + ^J {/c v / (o; 2 + 2/ 2 )}^; *+/>0. (32) 

OX J o 

This suggests at once suitable forms for the next approximation and for the 
free surface ; the equations are similar to (6) of the previous case, and we 
obtain in the same way 

^ = D-D 1 + X 1 , (33) 

where Di is a doublet at the image point (0, 0,/) and 



_ ™3 

dx 



Xi = ca rs 



\/(g/c)e- K ^- z hlfc e- 1 st xu Jo[/c^{(x + cu) 2 + y 2 }]8m(/cYu)du. 



(34) 
The corresponding surface elevation is 



•30 



7j = a? \ e- ,cf jQ{fc<\/(x 2 -\-y 2 )}fcdfc 



'CO 



'00 



+ a s \ \Z(gfc)e" K f 'tcd/c e~& u 3$[K^{{x + cu) 2 + y 2 }]%m(iiVii)du. (35) 
Jo Jo 

The first term represents the effect of the doublets D and Di. It can be 
verified by approximate methods that the second term includes a main part 
like the well-known wave pattern for ship waves. Since the expression in 
(35) gives finite and continuous values for the surface elevation, it might be 
of interest to examine some points in detail ; for instance, the elevation near 
the lines corresponding to the lines of cusps for a moving point disturbance. 
However, we pass now to the calculation of the resultant horizontal pressure 
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on the sphere. We have to find X 2 the image of Xi in the sphere ; for this 
we first put Xi into a different form by using 

7rJ [tf{(^ + ^) 2 -f 3/ 2 p] = cos {k (x -\- m) cos <f>} cos (Ky sin ^>)d<p. (36) 

Jo 

From (36) and (34), after carrying out the integration with respect to it, 
we obtain 



»oo 



ttXi = ca? \ e-«f-*) xdfc 









{ A sin (kx cos <£) + B cos (kx cos <p) } 

x cos (kij sin <£) cos <£ clef), (37) 



where A and B are given by (14) after writing c cos cf> for c. 

For convenience in the following analysis, we transfer the origin to the 
centre of the sphere, noting that in (37) we shall have exp. ( — 2k/+kz) in 
place of exp. ( — k/-{-kz). Also we use polar co-ordinates 

x = r cos a ; y = r sin a cos ft ; z — r sin a sin ft. 

The conditions for X 2 are that it must be a potential function, the 
disturbance due to it must ultimately vanish as we recede from the sphere, 
and on the sphere 

a(X 1 + X 3 )/3r = 0. (38) 

To avoid repetition of expressions like (37), we take out of it a typical term 

and write 

Xi = e KZ sin (kx cos (j>) cos (icy sin 0). (39) 

We know that the function 

T ~i e Ka*zir* s i n ( Ka 2 x cos ^J^) cos ( K0 2y s i n (^jf^ (40), 

satisfies the first two conditions for X 2 , but we find it does not fulfil (38).. 
An additional term is required, and it can be found in the following way.. 
Suppose that on the sphere we have 

e KZ sin (kx cos <fr) cos (icy sin <f>) = S A w Y m (a, ft), (41) 

where the right-hand side is an expansion in surface spherical harmonics.. 
Then for the term (39), all the conditions for X 2 are satisfied by 

ar -\p&zli* s i n ( Ka 2 x cos ^j r 2^ cos ( Ka 2y s i n ( j > /r 2 ) — %a m+1 A m Y m /(m-{- 1) r m+ K 

(42) 
Suppose, similarly, that on the sphere we have 

e KZ cos (kx cos (j>) cos (icy sin <£) = SB m Y wl (a, /3). (43) 

VOL..XCIII. — A. 2 S 
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Then the complete expression for X 2 is 



7rX 2 



f oq 



cw- 



e 2lc f icdtc 



m 7T 

ar -i e Ka?ziv* cos (xcPy s i n (^jr 2 ) cos </> 



x {A sin (km 2 x cos $>jr 2 ) -f B cos (/ca 2 x cos $/r 2 ) } d<j> 



— CCl' 



00 



JO 



f 75 " 

2 (AA W + BB ro )(m + l)" 1 (a/r) in+1 Y m cos <£ a 7 ^. 



o m 



We have now 



(44) 
(45) 



4> = D-Di + Xi + Xa = D+X, 
and the pressure equation is 

jj/p = —cd<f>/dx—gz + /j,(f> — ^q 2 . (46) 

The wave resistance, or the resultant horizontal pressure on the sphere, is 



E 



"IT 

d 



a 



•2 



a 2 p sin a cos a d/3. 



m 



Omitting terms which, from symmetry, will give no contribution to B, we 
have 

i 3d ax 



p dx dr dr r 2 doc doc r 2 sin 2 a 9/3 9/3 



(43) 



But when, r = a, we have 

3D/9^ = ; 9D/9a = -\m sin a ; 9X/9r = 0, 
hence p/p = (3c/ 2a) sin ocdX/doc-}- fiX. (49) 

We must, now substitute (49) in (47) and use the value of X given by the 
sum of (37) and (44) on the sphere ; it is clear that we may omit the doublet 
Di as it will not affect the limiting value of B when p, is zero. 

6. Consider, in the first place, the contribution of the first term in the 
value of p given in (49). In the repeated integrals which are obtained, we 
may change the order of integration, and we shall carry out first the summa- 
tion over the surface of the sphere. We notice that, when r = a, the first 
term in the value of X2 in (44) is equal to the- value of Xi ; the additional 
part of X2 is the term involving the expansions in spherical surface 
harmonics.. Choose a typical term from the latter part, and we find we have 

to evaluate 

J|sinacosa(9Y m /9a)^S, (50) 

taken over the surface of unit sphere. 

But this integral is equal to 

-3JP 2 (cos«)Y m («, /3)dS. (51) 

Hence, the only term which has a value different from zero is the term in 
Y 2 , the surface harmonic of the second order. From the manner in which 
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the expansions were introduced, in (41) and (43), it follows that the contri- 
bution of the second term in (44) is one-third of that of the first term ; 
hence, summing up the result so far as the first term of (49) is concerned, we 
have 



00 







IT 



C2it rir 

ttE' = -5cVp e- 2K ffcdfc co3<£d<£ dfi\ sin aP 2 (cos a )e Kashlasin ^ 

o 



o Jo 
x cos {tea sin a cos /3 sin <f>) 

x {A sin(/cct cos a cos </>) + B cos (/eacos a cosc£)} da. (52) 
Taking the integration with respect to /3, we find it is equal to 



IT 



2 gicasmacosp cos ^ Ka s { n a s j n ^ gm ^ <^g _- 2tt Io {jca sin a cos <£), (53) 

Jo 

where I (x) is the Bessel function J (ix), a result which may be obtained by 
direct expansion and integration term by term. For the integration with 
respect to a the term in A in (52) obviously gives zero, and we are left with 



'IT 



2w Io {tea cos (f> sin a) cos {tea cos </> cos a) P 2 (cos a) sin a da. (54) 

Jo 

Here also we may expand in powers of act and integrate term by term ; 
it can be shown that the integral of the coefficient of (tea) 71 vanishes except 
for the single term tc 2 a 2 ; thus we find that (54) reduces to 

— (2tt/ 5 ) k 2 <z 2 cos 2 (£. 

7. We have now to consider the term fiK in the value for p in (49). We 
might omit this term, on general grounds, as giving no contribution to E 
ultimately when fi vanishes ; for X is the velocity potential for a sphere at 
rest in a given field Xi. However, it may be left in, and we have a similar 
calculation. Taking the second integral in (44), we find it is now only the 
term in Yi which counts ; hence the contribution of this part is one-half of 
that of the first integral in (44). Further, it is the term involving A which 
gives a value different from zero when integrating with respect to a, and 
instead of (54) we have 



fir 

2tt 
Jo 



I (/co, cos cf) sin a) sin (tea cos $ cos a) P x (cos a) sin a da, 



which reduces to (47r/3) /ca cos $. 

8. Collecting the various results, we have now 

E = — 2ca 6 p\ e- 2K f K 2 d/<: (kcBcos <j> + pA) cos 2 <j>d(b, (55) 

Jo Jo 

a form which may be compared with the corresponding expression for the 
cylinder in (29). 
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A and B are given by (14) when we replace c by ccos<£; putting these 
values in (55), we see that we may change the order of integration. Further, 
as we make fi vanish ultimately, we may use simplified forms of A and B 
corresponding to (16). These give 

Jo Jo (^—^osec 2 ^) 2 + (/x/c) 2 sec 2 

To obtain the limiting value for fi zero we may treat this like the similar 
expressions in (30); or, alternatively, we may put (/x/c)sec<£ = 1/n, and use 
the general result 



Lim 

n — > <x> 



,l+n 2 (x — a) 2 2 u v /'t/v n 



The apparent difficulty with regard to values of <fc near tt/2 is overcome by 
noticing that with the particular functions involved in II no extra contribu- 
tion arises from such terms near the upper limits of the variables. Carrying 
out the integration in k in this way, and changing the remaining variable by 
putting tan <fi = t, we obtain 



•00 



E = 4:irc/pa?c-h- 2 ^ \ (l + t 2 y; 2 e~ 2 ^ c \ (56) 

The remaining integral can be expressed in terms of known functions. 
Possibly the simplest method is to use the confluent hypergeometric 
function* defined, for real positive values of a and for real values of k and m 
for which k — m — -| <^ 0, by 

e ~a{2 k fee 

W k>m (*) = - r ( 1 _ h + m) j q u-*-* +m (1 + u/*y>-i+ m e-» du. (57) 

We have now the wave resistance of the sphere given by 

K = \irWgpaSf-Zo?l 2 e-*l 2 W U1 (a); a = 2gf/c 2 . (58) 

8. For purposes of calculation, we require expansions of Wi,i(a). This 
function belongs to the logarithmic type of confluent hypergeometric function, 
and general expansions are not available in this case ; however, they can be 
obtained without difficulty for Wi,i. In the first place, the differential 
equation satisfied by Wi,i is 

d ^ + (-i + l-A\ y = o. (59) 



da 2 \ 4 ol 4 



We use the ordinary methods for solving by means of power series. The 
roots of the indicial equation are -f and — J ; hence one of the fundamental 

* E. T. Whittaker, 'Bull. Amer. Math, Soc.,' vol. 10, p. 125 ; also Whittaker and 
Watson, ' Modern Analysis, 5 Chap. XVI. 
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solutions will contain logarithms. Calculating the coefficients step by step ? 
we obtain as a fundamental system 

We know that Wi,i is a linear function of yj and 3/2; however, it is simpler to 
obtain an expansion directly and use (60) to verify it. For this purpose we 
use the equivalent contour integral for the confluent hypergeometric function, 

r(s)r(-s-^m + J)r(-s-fc + m + J) 



* » /(,•> to ■ — 



2tt£ 



,00* 



F(~ &— m + J)r(— k + m-\-l 2 ) 

where the contour has loops if necessary, so that the poles of T(s) and those 
of r(— s— A— m + J)r(— s— Zs + ^H-J) are on opposite sides of it. The 
integral can be evaluated by the method of residues. When k = m = 1, the 
poles at which the residues have to be found are simple poles at s = — |-, — f , 
together with double poles at 6' = -|, -|, -§-, .... The latter series gives rise to 
logarithmic residues. Carrying out the calculation, we obtain 

Wl;1 = w -^-w ( a -, 2+ | ._„) _ 4 1 A - a/2 | log a I T ^±i^ * 

+ ^( 7 -21 g 2 -g) + S^| r(p r(pH-g + 8 J > (62) 

where 7 is Euler's constant 0*5772.... The coefficients may be put into 
alternative forms more suited for calculation : for instance 

dpT(p + l)T(p + Z) 

2? .p ! (jp -f 2) ! L to T n (2^ — 1) 7 ?i 

For numerical calculation we have 

-(^- 1 °«i a )( 1 + r + 3 L 2 a2 + 4 a3+ -)- (63 > 

The expansion may be confirmed by comparison with the fundamental 
solutions of the differential equation given in (60) ; we find that 

(8/3)7r*W lfl - {2\og2-v~^ yi ~y 2 . 

For large values of a the general asymptotic expansion of W k , m is available ; 
and in this case we have 

Wi,i~««-^fl + |- + ^^-# a -,+ -)• (64) 

4 a M of 128 a* I v 7 
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532 Some Cases of Wave Motion due to a Submerged Obstacle. 

9. With (63) and (64) we can now calculate the resistance E from (58). 
For a given depth/, the variation of the resistance with the velocity is shown 
in the following curve, for which E has been calculated for various values of 

clx/(gf\ 




j 



The curve is very similar in form to the two-dimensional case of a circular 
cylinder. For small velocities, that is a large, if we take the first term of the 
asymptotic expansion (64), we have 

E = x /(27rY/f) • pcPc-*e-W* 9 

which may be compared with (30) for the cylinder. It is of interest to 
notice the similar law of variation of wave resistance with speed for the few 
cases of rigid bodies which have been worked out. The method adopted here 
can be applied to bodies of different forms, and it is hoped to illustrate later 
some interference effects. 
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